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The problem of the equihbrium state of the charged many-particle system above dielectric surface is formulated.We 
consider the case of the presence of the external attractive pressing field and the case of its absence. The equilibrium 
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1 Introduction 

The problems concerned with the research of the 
charges above dielectric surface belong to classical elec- 
trodynamics and electrostatics problems. A special 
interest to such problems appeared due to the phe- 
nomenon of the Wigner crystallization. These re- 
searches were initiated in 1934 by Wigner in his theo- 
retical work [T], per se. In this work, the possibility of 
the existence of periodic structures in the systems with 
repulsive forces between particles was demonstrated 
by the example of the crystallization of the three- 
dimensional low-density gas of electrons in the field 
generated by the spatial-homogeneous positive charge. 
This field played exactly a role of the compensative fac- 
tor for the repulsive forces. The Landau-Silin Fermi- 
liquid theory also enables to predict the existence of 
spatially-periodic state of electrons in metals and to 
describe its structure (see in this case Ref. [5]). The ex- 
perimental improvement of Wigners' prediction of the 
three-dimensional crystallic structures still does not ex- 
ist (see, e.g., Refs. [31 II]). This is caused by difhcul- 
ties in the achievement of the experimental conditions 
for the mentioned phenomenon, which is also called as 
"Wigner crystallization". 

However, as it is well known, different two-dimension 
periodic electron structures above the surface of a 
fluid helium are experimentally realized (the so-called 
" Wigner crystals" ) . The works O El El [8] may be re- 
ferred as the first publications containing theoretical 
and real experimental results of different properties of 
the surface electrons. The large number of works re- 
lated to the theoretical and experimental research in 
this area has appeared by now. 



The theoretical papers that are devoted to the mi- 
croscopic description of the charge state above dielec- 
tric surface are usually based on the conception of 
an isolated charge above dielectric surface interact- 
ing with its electrostatic reflection in dielectric ("lev- 
itate electron", see, e.g., Refs. [3l IH |9]). In this 
case, the quantum- mechanical state "charge - electro- 
static reflection" is described as the hydrogen like one- 
dimensional state with the corresponding energy struc- 
ture. Very often the localization of such quantum- 
mechanical object in the ground state is considered 
(see Refs [3l|4l|5l|6l|9]) occuring at some distance from 
the surface (first "Bohr radius"). This, particularly in 
most cases, allows not to take into account the influ- 
ence of the surface inhomogeneity on the single charge 
state. However, at the description of many-particle 
charge system above dielectric surface the mentioned 
approach inevitably faces some difficulties. For exam- 
ple, such difficulty appears when the electron density 
above dielectric surface does not allow to consider the 
charged particles as isolated, i.e., it is necessary to take 
into account the interparticle interaction. 

The references, which are devoted to the two- 
dimensional Wigner crystallization in the phenomeno- 
logical approach, predominantly consider the system 
that consists of a large number of charged particles near 
the surface of the fluid dielectric as a two-dimensional 
structure (see, e.g., Refs. [3 HEl 012118]). 

Basing on the premises, it becomes clear that the 
complete description of charges above dielectric sur- 
face needs to take into account their spatial distribu- 
tion in vacuum. The possibility of charges adsorption 
by the surface must also be taken into consideration 
(in this case the surface inhomogeneities must play a 



crucial role). The possibility of charge spatial distribu- 
tion above dielectric surface comes from the fact that 
a charged particle is always attracted by a dielectric 
surface. Moreover, in the experiments [31 Uj concerned 
with the registration of two-dimensional Wigner crys- 
tallization an external electric field attracts charges to 
the surface and affects on their spatial distribution. 

The present paper is devoted to the problem of equi- 
librium charge distribution above dielectric surface as 
in the external pressing electric field as in its absence. 
This problem is considered in the case of ideally plane 
vacuum-dielectric boundary and in the case of "wavy" 
spatial-periodic surface with account of the possibility 
of existence of the"sticked" to the surface charges. In 
our opinion, the formulated problem is interesting as 
from purely academic side as from the research side 
of the influence of volume charges located closely to 
the fluid helium surface on the spatial-inhomogeneous 
state of the charges adsorbed on the hclimn surface. 



2 Equations of electrostatics for 
many- charge system above di- 
electric surface 

Let us consider the equilibrium system of charged par- 
ticles (Fermi-particles) with the charge Q per particle 
that is situated in vacuum above dielectric surface with 
the permittivity e. We describe below the surface pro- 
file by function ^(p) = ^(x, ?/), where p = {a;, y\ is the 
radius-vector in the plane 2 = of Cartesian coordi- 
nates {x, y, z}. The vacuum - dielectric boundary lies 
in the plane z — and we consider it unbounded below. 
All physical quantities considered in the area above di- 
electric, i.e., at z > ^(p) we mark by the index " 1" and 
all quantities concerned with the dielectric {z < £,{p)) 
we mark by the index "2". Let us assume that the 
external pressing electric field E acts on particles and 
is directed along the z-axis. We also assume the exis- 
tence of some potential barrier that forbids the charges 
to penetrate inside the dielectric. 

As it is mentioned above, the charged particles are 
always attracted by the dielectric. Therefore, even in 
the absence of the external pressing electric field there 
is a reason to believe that there are conditions under 
what the stable equilibrium distribution along z-axis 
is developed. To avoid the questions on the repulsion 
of likely charged particles along the plane p we shall 
consider the system located in a vessel with the walls 
at p ^ 00. These walls forbid the charges to leave the 
system. Let us describe the equilibrium charge distri- 
bution above the dielectric surface by the distribution 
function /(p; z, p). 

The electric field potential ipi in vacuum above the 
dielectric surface must satisfy the Poisson's equation 



where A is the Laplace operator, 

+ ^""a^^ + a?' (2) 

9{z — C(p)) is the Heaviside function. In eq. ^ the 
quantity n(z^ p) is the charge density above the dielec- 
tric surface, which can be expressed in terms of the 
distribution hmction /(p; z, p) as 



n{z, P)^ j d^PfiP', p) 



(3) 



As charges are considered as Fermi-particles, the dis- 
tribution function /(p; z, p) has the following form: 



/(p;^,p) 



X < exp (3 



{2nhf 

,2 



^-f Q(/3l(z,p)-^ 



(4) 



where g = {2Sq -1-1), Sq is the spin of the charged 
particle, f3 = T is the temperature in the en- 

ergy units, m is the charge mass and /i is the chemical 
potential of charges. We emphasize that taking into 
account relations ([3]) and (j4|) the eq. ([T]) is often called 
the Thomas-Fermi equation. 

The electric field potential (p2 is the result of charge 
absence in the dielectric. In the assumption of the 
dielectric homogeneity and isotropy it must satisfy the 
Laplace's equation 



eA(p2{z,p) = 0. 



(5) 



If the system is placed in the external static homo- 
geneous electric field, the potentials (fi and ip2 can be 
written in the form 

Vl = y^l^ + ^[""^ , V2 = + "Jf : (6) 

where '~p\' ^ ^^2 potentials induced by the sys- 

tem of charges in vacuum and in the dielectric, re- 
spectively, Lp\'' and ip^2^ are the potentials of the ex- 
ternal field in vacuum and dielectric. According to 
eqs. ([I]), (O, these fields satisfy the following equations: 



A(^« (z, p) = -47rQn(z, p), A^^"^^ = 0, 
A^W=0, A^(^)=0. 



(7) 



A(pi(z, p) = -47rQn(z, p)d{z ~ ^(p)). 



(1) 



Eqs. (dl), ([5]) for the potentials must be expanded 
with the boundary conditions on the vacuum-dielectric 
border (as usual, these conditions can be obtained di- 
rectly from eqs. ([I]), dS]), see e.g.TOj): 

¥'i(^;,p)U=^(p) = ¥'2(2;,p)U=4(p), 

ni{p){e^^Lp2{z, p) - £Vj(^i(z, p)}2=^(p) (8) 
-47ra(p,e(p')), 

where n(p) is the unit vector of the surface normal in 
the point p, ct(p, ^(p')) is the surface charge density 



in the point p (here we emphasize the functional de- 
pendence of this value on the surface profile ^(p') ). 
The surface charge density a{p,^{p' j) must satisfy the 
following relation: 



J dS^a{p,ap'))^QN^, 



(9) 



where A'^ is the complete charge number on the dielec- 
tric surface and dS^ is the surface element with the 
profile £,{p'): 



(10) 



The surface charge density can appear due to several 
reasons. For example, these charges may be specially 
placed on the dielectric surface and can stay there for 
arbitrary long time. At that time, the surface charges 
and charges above the dielectric surface can differ in 
sign. But in this case we need to consider the possi- 
bility of the formation of bound states of the opposite 
charged particles. Taking into account the presence of 
such bound states represents a separate rather com- 
plicated problem. In the present paper such case of 
surface charges is not considered. The case, in which 
some part of charges condenses on the surface from the 
volume distribution, stays for some period of time and 
then back to the volume, is possible too. In this case, 
the equilibrium distribution of charges above the sur- 
face that coexist with the " sticked" surface charges for 
some period of time (the lifetime of the charge staying 
on the surface) is possible. Next, we shall take into ac- 
count only the possibility of the presence of the surface 
charges that have the same sign as the volume ones. 

It is easy to see that in eq. ^ the directional cosines 
of the surface normal vector n(p) in the point p play 
the main role: cos i/ is cosine of the angle between the 
normal and z-axis, cos A is cosine of the angle between 
the normal and x-axis and cos n is cosine of the angle 
between the normal and y-axis. In the case when the 
surface profile is given explicitly (in our case z = C(p))j 
these cosines are determined by the following relations: 



cosi^ = 



cos A = 



cos fj, ■ 



1 




Hdap)/dpr' 




dap)/dx 




+ mp)/dp)^ 




di{p)/dy 




+ {dS,{p)/dpf 



(11) 



The derived electrostatics equations ([T]), ([5]) with the 
boundary conditions ([8]), pTjl can be solved analyti- 
cally in a very low case count. Some of these cases 
are considered below. Before solving eqs. H]), ([5]), let 
us consider the simplification of the boundary condi- 
tions ^ in the case, when the surface profile differs 
a little from the plane one. In this case we essen- 
tially have the effective boundary conditions. From 



eqs. ([TUl) . (fTT|) it is obvious that the surface slightly dif- 
fers from the plane one, when the surface profile slowly 
varies on coordinate, i.e., when the following inequali- 
ties take place: 



mp)/dx\ « 1, mp)/dv\ « 1. 



(12) 



Let us also consider that the surface profile ^(p) can 
be given as: 



e(p) = e + e(p), l(p) 



It is easy to see that the inequality H 
this case by the conditions 



(13) 
is provided in 



dl{p)/dx « 1, dl{p)/dy « 1 



(14) 



The directional cosines ([TT|) with accuracy up to the 
second order over 9^ (p) /dp have the following form: 



cosvK,!, cosX — ~d({p)/dx, 
cos/i -dl{p)ldy. 



(15) 



If the relations ([T^l - P^ take place, we can expect 
that the charge and the field distributions in the system 
slightly differ from the distributions that take place 
in the case of the plane dielectric surface. Then, the 
potentials ipi{z,p) and ip2{z,p) (see eqs. ([1]), (O) can 
be written as 



<fii{z,p) ^ ifii{z) + ifii{z,p), 

^2(2, p) ip2{z) + ^2(2, p), 



(16) 



where ^pi{z) and ^P2{z) are the potentials of some elec- 
tric field above the dielectric and inside of it (but not 
on the surface!) in the case of the plane surface. The 
small distortions of the field above the dielectric and 
inside of it are described by the potentials Lpi {z, p) and 
(p2{z, p) due to the surface inhomogeneity in the men- 
tioned above sense. The meaning of the introduced po- 
tentials (/ji(2;) and ip2{z), and also ipi{z,p) and ip2{z, p) 
becomes more clear after the obtaining of the Poisson's 
equations and effective boundary conditions for them. 

According to the assumption of small field pertru- 
bations provided by the wave surface, the following in- 
equalities take place: 

l<^l(^)l > l<^l(^,P)h l'^2(^)| > |<^2(2,P)|. (17) 

Let us also consider that the distribution of charges 
that can be condensed on surface slightly differs from 
the homogeneous one: 



cT{p,0^^{0 + ^ip,0 + 

1^(01 »|a(p;OI, 



e(p), 



(18) 



In the expressions (|18p the quantity ct(p, ^) corresponds 
to the impact of the weakly inhomogeneous charge dis- 
tribution on the plane dielectric surface with z — ^ pro- 
file. The quantity ^{p)[dcr{£,)/d£,] in eq. ([18]) describes 



the surface charge inhomogeneity related to the weak 
irregularity of the surface itself. 

In the expressions (|16p . (fT7)) we consider that in 
the case of uniformly charged surface, which is ideally 
plane and infinitely extended with charge density cr(^) 
both the field and the charge distributions are homoge- 
neous along p plane. In other words, the spatial charge 
and field distribution depends only on z coordinates. 

From the relations (fT51) - (fT5)) it is easy to obtain the 
effective boundary conditions for field potentials on 
vacuum-dielectric boundary in the case, when the di- 
electric surface weakly differs from the ideally plane. 
To this end, we must develop the perturbation theory 
over small values f(p), and d^{p)/dp, which, 

according to the expressions (fT2)) - lfT8|) . can be given as 
follows: 



d ^ d ^ 

e — [ip2(.z) + ip2{z,p)] - -^VPiiz) + ^i{z,p)] 
oz oz 



= 47ra(p;^ + e(p))- 



(19) 



Making the necessary calculations up to the first order 
of the perturbation theory from the first relation of 
eq. we obtain 



Lpi{z%=(^ = </72(2)U={, 
^l{z,p) - ^2{Z,P), 



(20) 



9(^2 (z) ^^pl{z) 



dz 



dz 



The use of the perturbation theory up to the first or- 
der for the second relation of eq. results in the 
following equalities: 



dz 



dz 



2=? 



Let us remind that in electrostatics the denotations like 
{d'^(pi{z)/dz^}z=^, {d-^ip2{z) /dz^} have the mean- 
ing of limits 

{d^^i{z)/dz^},=^ = lim{aVi(^)/9z'}.=?+/., 

{d^Mz)/dz^},=^ - lim{9V2(2)/5z2}.=C-h. 

Then, for the further simplification of the obtained ef- 
fective boundary conditions ([20|) . ([2T|l . according to 
eqs. dl]), (ID), (fTBjl we can use the following equations, 
which are satisfied by the potentials ipi{z), (p2{z): 

|J^i(z) = -4^Qn(z)0(z-O, ^Mz) = 0, (22) 



where 



f{p;z) = 



n{z) = J dV(P, z) 



(27rft)3 



exp f3 



- — h Qipi{z) - fi 

2m 



(23) 

Taking into account eq. pT|) , the conditions (pO|) can be 
expressed in the following form (the first one of them 
remains the same): 



_dip2iz) d(pi{z) 



dz 

Att-I Qn{z) - 



dz 
da{0 



di 

d^i(z,p) 
dz 



= 47ra(0, 

i{p)-A^a{p-0= (24) 
^d^2{z,p) 1 



dz 



Thus, we obtain the effective boundary condi- 
tions (PO)) . (|24p for the fields in the system of charges 
above the dielectric surface with the surface profile that 
slightly differs from plane 



+ Apipi[z,p) = 47rQ 



d'^^2iz,p) 
dz"^ 



d^ 

Apip2{z,p) = 0. 



'Pi{z,p), 



(25) 



3 System of charges above the 
ideally plane dielectric surface 

It is easy to see that the obtained equations ([^ . (^5)) 
of electrostatics and the effective boundary condi- 
tions (IM)) for these equations are much simpler 
than the initial electrostatic equations ([T|) , (O and the 
boundary conditions ([8]). Firstly, to solve the equa- 
tions that determine the charge and the field distribu- 
tion above the vacuum-dielectric boundary one needs 
consider the case of ideally plane surface of this border 
that lies in z = ^. Let us start withe considering the 
case of the surface in the absence of the charge cr = 0. 
Then, the solution of (|22p must satisfy the following 
boundary conditions: 



(pi(z)U=^ = (P2iz) 

_dip2iz) d(p2{z)) 



dz 



dz 



0. 



(26) 



To simplify the further calculations, let us write the 
first formula from eq. in the following form: 



d^liz) 



oo 

= -47rQi/ J dee^/^{ exp/3(e - i/') + 1}, (27) 



where we denote 

V^(z) = ^-Q(pi(z), iy= {2m)^/^/2TT^h^. (28) 



Here we also consider the spin of a charged particle 
equal to 1/2, ip is so-called electrochemical potential. 

Multiplying eq. ^7} by the derivative {d(pi{z)/dz) 
and using the following equality 



In 



1 



after simple calculations we obtain the first-order dif- 
ferential equation: 



2 °^ 



where C is an arbitrary integration constant. Thus, 
the need of the following equation solving arises: 



dz 







1/2 



(29) 

The sign before the square root in eq. (^5)) must be cho- 
sen from the following consideration. The force acting 
on the charges at 2: > ^ presses these charges to the 
dielectric surface. Thus, in the case of positive charges 
above the dielectric we choose the positive sign, and 
in the case of negative charges we choose the negative 
one. Let us consider below the distribution of negative 
charges above the dielectric surface, Q = — e, e > 0. 
Hence, the potential ipi satisfies the relation: 



00 


Now we make the following denotations: 

ipi{z = 0) = ipo, i>{z = 0) = ^ + eipi{z ^ 0), 
f d^i{z) \ 
V dz 



1/2 



(30) 



(31) 



Let us remind that we consider the case of the elec- 
tric forces that attract charges to the dielectric surface. 
Thus, at z ^ CX3 there is no charges, /(p; z) — > 0, or 

Z — 'OO 

{exp/3(e- V) + l}"^ ^ 0. (32) 

Z — >OQ 

The action of the electrostatic image force along z-axis 
must vanish at z 00: 



dz 



0. 



As the result, it is essential to say that at z — + c» the 
following relation takes place: 



dz 



a = 

oz 



(33) 



where E is the external field intensity that attracts 
charges to the dielectric surface. 



At z = from eq. ([50)1 one can get 

00 

El = j de£3/2 {exp/3 (e - ^q) + 1}-^ + C. 


On the other hand, from the same equation and taking 
into account eqs. (32), (33) at z — > 00 we obtain: 



C = E^ 



(34) 



Comparing the last two expressions, we come to the 
relation between the constants -00 1 Eq (see eq. (I5T]) ') 
and the external electromagnetic field E: 

00 

El~E^ = j de£3/2 {exp/3(e-7Ao) + l}"'. 


(35) 

Then, after integration of the first expression from 
eq. (|23p over z within the limits from ^ to e and us- 
ing eqs. (1211), dSll), we get: 



Eq 



Aliens, e > 0, 



(36) 



where Ug is the number of the volume charges per unit 
of the plane dielectric surface: 



Us = J dzn{z), 



00 

n{z) = i' J dee^/^ {exp (3 {e - ip{z)) + iy\ 




(37) 



Let us emphasize that for the equilibrium charge sys- 
tem above dielectric the value of the number de- 
pends neither on the coordinates, neither on the fields' 
distribution. It is determined only by the entire num- 
ber N of the charges above dielectric. We also point 
out that this value characterizes the additional field 
intensity that presses the charges to the dielectric sur- 
face. Besides that, this field is generated by the charges 
themselves. 

Thus, eqs. ([55)) . ([M)) allow to express the unknown 
quantities ipo and Eq (integration constants of eq. (|27p ) 
in terms of the external pressing electric field E and the 
number of charges above the unit item of the dielectric 
surface rig (see eq. (|57)) ). 

The second equation in (|38p can be solved trivially 
in general case, because the electric field intensity in 
dielectric does not depend on z. Using the boundary 
conditions (j26p . we can express the potential of electric 
field in dielectric in the following form: 



Eo 

tp2 = Z ■ 

e 



(38) 



where E2 is the electric field intensity in the dielectric, 
Eq can be expressed from eqs. (|35t . ([36)) and (po is the 
potential on the surface. As in the electrostatic case. 



the field potential is determined accurate within a con- 
stant. Therefore, we set the potential (po equal to zero 
below. In this case, the value of the electrochemical 
potential on the dielectric surface coincides with the 
chemical one: 



Ipiz = 0) = IpQ = ipQ = 0. 



(39) 



Taking into account eq. ([M)) . the spatial distribution 
of the potential (see eq. (|30p) can be written as follows: 



dtpi 
dz 



167r 



x{exp/3(e- V') + 1} ^ + E^' 
ip{z) EE ^ + eipi{z). 



1/2 



(40) 



It is easy to see that in general case the solution of this 
equation can be found only in quadratures (see below). 
However, the gas of charged Fcrmi-particles above the 
dielectric surface is nondegenerate. Therefore, the so- 
lution of eq. (|40)l can be obtained analytically. Indeed, 
in the case of nondegenerate gas its distribution func- 
tion has the form that weakly differs from Boltzmans' 



{exp (3 {e — ip) + 1} ^ ~ exp (3 {ip — e) . 

Accordingly, the expression for the density distribution 
of a gas along the z coordinate (see eq. ([37])) becomes: 

n{z) « ^i//3-3/2 exp(/3V). (41) 

As the Fermi-particle gas is degenerate at low tem- 
perature and high density ranges (see e.g. [11]) from 
eq. (|4ip one can get the gas nondegeneracy condition: 

exp(/3V) < 1, 

As the electrochemical potential depends on z, this con- 
dition is obviously realized in the case when the follow- 
ing inequality takes place: 



expiPipo) < 1, 



(42) 



where "00 is the electrochemical potential on the dielec- 
tric surface (see eqs. (pS)) . ([M]) in this case). The last 
statement takes place because of the assumption of the 
particle absence at z ^ oo, see above. So, according 
to eq. ([4T|) . ([42]) . the formula ([40]) can be expressed as: 



1/2 



|j = - {4^3/2^2^-5/2^ g^p(^^) ^ g2^2| 

This equation has the analytical solution: 

f.r^/^exp(/.^(.)) = 4^^, 



(43) 



(44) 



^Piz) 



eipi{z), 



where the function xi^) is defined by the relation: 

Eq — E 



X{z) 



Eo + E 
zo = {l3eE) 



exp{-(z-^)/zo}. 



(45) 



Let us emphasize that the multiplier before the expo- 
nent in eq. P5]l according to eq. (|36p can be expressed 
in terms of the intensity of the external electric field 
and the number of charges rig in the "column" above 
the surface unit element: 



Eq-E 
Eo + E 



E + 2Trens 



From the eqs. ([4T]), ([43]) and ([44]) follows that the 
charge density above the dielectric surface has the dis- 
tribution: 

.E^ ix{z) 



n{z) = (3- 



(46) 



and the electric field intensity above the dielectric 
Ei{z) is expressed as 



Ei{z)^E 



i + x(^) 



(47) 



It is easy to see that at high values of z 3> (see 
eq. ([45]) ). the charge distribution above the dielectric 
surface is close to the Boltzman distribution and the 
electric field density exponentially tends to the external 
pressing electric field density. This fact confirms the 
above assumptions (see eqs. ([5^ . ([33])). 

The inequality ([1^ that determines the nondegener- 
acy condition of the charge gas can be written in terms 
of the obtained solutions: 



ensv~'^(3'''/'^ {E + 2iiens) < 1. 



(48) 



It is obvious that this inequality is not accomplished 
in the case of low temperature range or high values of 
the external pressing field. 

Expressions ([^^ - (|15)) allow to make the limit pro- 
cess at i? ^ 0. In the case of the absence of the external 
pressing field, these solutions have the following form: 



Ei{z) ^ £^0 



1 



p2 r 



z-^ 
2zo 

-2 



2zo 



Ea/e, 



where 



Zo = {PeEoY 



Eq = 'ineus 



(49) 



(50) 



If to compare the expressions ([46]) . ([47]) and ([49]) . it 
is easy to see that in the case of the absence of the 
external pressing field the exponential law of the elec- 
tric field and charge density above the dielectric surface 
changes to the weaker power depedence. In this case, 
the inequality ([^5]) can be written as: 



(51) 



(en,)'i/-i/35/2«l. 



Note, that it takes place in the region of the relatively 
high temperatures and low charge number in the vol- 
ume above a surface area unit, see eq. p7p . 

The obtained formulae are the so- 

lution of the problem of the field and nondegenerate 
charged gas distribution in charged particle system 
above the plane dielectric surface as in the external 
pressing field as in its absence. Let us emphasize that 
the dielectric permittivity does not appear in these ex- 
pressions. The reason is that the problem is homo- 
geneous along the surface coordinate p. In the case of 
inhomogeneity along p, the solution of equations essen- 
tially depends on the sort of the dielectric, i.e., on its 
permittivity e. These inhomogeneities may be caused 
by the inhomogeneities of the surface itself or by inho- 
mogeneity charge distribution on it (or the both rea- 
sons simultaneously, see eqs. ©, (HH), (UHl), (EH))- 

In the case of the degenerate gas, i.e., when the con- 
dition (HS)) or (|?T|) fails, the solution that is obtained 
earlier is inapplicable. Let us make the following re- 
mark relating to this fact. As it is mentioned earlier, 
the charge density distribution decreases with the dis- 
tance from the surface. For this reason, in general case 
described by eq. (I40p the gas can be degenerate in the 
area near the dielectric surface and nondegenerate far 
from it. The typical distance from the surface that sep- 
arates these cases can be obtained using the following 
considerations. As it is well known (see e. g. Ref. 
in low temperatures region the temperature expansions 
are widely used for the calculus of the thermodynam- 
ical quantities characterizing the gas. Applying such 
expansion to the integral over the energy in eq. ()40|) . 
we obtain: 



oo 



-1 



(52) 

From this expression it is easy to see that such expan- 
sion is absolutely useless near the point zi obtained 
from the condition 



ip{zi) = + eipi{zi) = 0. 



(53) 



The solution of eq. (PO)) obtained in quadratures is 
given by 



Z-^= -Zq 



pil^ oo 
0^0 



x{exp(2/-C) + l} +1 





-1/2 



(54) 



where the distance zq is determined by eq. (|45p . Tak- 
ing into account eqs. ([53)1 . ((54)) . the expression of the 



border distance zi can be written as 

I (55) 



:{exp (y-O + lj-' + lj 



-1/2 



where the electrochemical potential V'o as the function 
of temperature T = j3~^ and the external electric field 
is obtained from the equation (see eqs. ([35|l. ([36]), ([39]) ) 



{Airens + Ey - E^ = 



167r 



OO 

j dee^/^ {exp/3(e- i/>o) + 1}"^ 



(56) 



As it is mentioned above, the potential i^iiz) is defined 
accurate within an arbitrary constant, which can be set 
equal to zero. Hence, eq. (j56p is the expression defining 
the chemical potential /i, see eq. 

As expected, the typical distance zi (see eq. (155)) ') is 
defined by the temperature, the external pressing field 
and the number of charges above the dielectric surface 
area unit. Thus, the charge gas is nondegenerate in the 
region z 3> 2i and degenerate at z <C ^i- Let us point 
out that the solutions (|15|) - ([5(I)) are obtained assuming 
the charge gas nondegeneracy in the entire area above 
the dielectric surface. Therefore, in general case the 
mentioned expressions describe the charge system only 
in the region z 3> 2^1 . The charge gas above the dielec- 
tric surface can be degenerate even in the case of the 
absence of the external pressing field. It is easy to see 
if to analyze the expression ([55]) at £^ ^ with the 
account of eq. . 

In the case of the generate charge gas above dielectric 
surface, eq. (PD|) according to eq. ([5^ can be written 
in more simple form: 



dipi ( 32tt 
~d7 ^ " lis" 

Tp{z) = ^ + eipi{z). 



1/2 



(57) 



However, in this case eq. ([57]) can not be solved an- 
alytically, and the numerical integration methods are 
needed. 

Let us show now the influence of the plane dielec- 
tric surface charges on the obtained results in the 
present section of this paper. It is well known that in- 
finitely thin homogeneously charged plate with charge 
density cr(^) induces the homogeneous field intensity 
E^ = 2'Ka{^) in vacuum (in this case the expression 
(t(^) shows that the surface plane is described by the 
equation z = This field has the opposite direction 
in the opposite sides of the plane. In the case of the 
charged plane dielectric surface the situation is abso- 
lutely similar. E.g., the negatively charged dielectric 
surface induces the field intensity £^20- — ^tt |ct(^)| /e in 



the dielectric and E^^ — —2TT\a{^)\ above the dielec- 
tric surface (see eq. ([Ml))- As mentioned earlier in the 
present paper, we consider only the cases of the same 
signs of charges as on the dielectric surface, as in the 
volume above it (in our case we consider the negative 
charges). In this case, the field induced by the surface 
charges repulses the volume charges from the surface. 
So, the results obtained in the present partition remain 
useful if we substitute the external electric field in vac- 
uum E for E — 2Tr\a{£_)\ in the expressions ((3T|) - ((56)) . 

£; ^ 27r|cr(^)| . (58) 

It is easy to see that it is necessary to satisfy the con- 
dition 

Eo-27:\a{0\>0 (59) 

that provides the possibility of existence of the equilib- 
rium volume charge distribution above the dielectric 
surface in the repulsive field of the surface charges. 

4 The charge system above the 
spatially inhomogeneous di- 
electric surface 

As already mentioned, the spatial inhomogeneities can 
be caused by the surface heterogeneities or by the inho- 
mogeneous charge distribution on it (or the both rea- 
sons simuhaneously, see eqs. (0), (HH), ([201), (EH)))- Let 
us consider the mentioned surface inhomogeneities that 
slightly distort the electric field induced by the charge 
system above the plane dielectric: 

(pi{z,p) = ipi{z) + (pi{z,p), \ipi{z)\ > |(pi(z,p)| , 

<f2{z,p) = (P2{z) + <f2{z,p), \V2{z)\ > \(f2{z,p)\ , 

where 1^31(2:), 952 (-z) sue the potentials above the dielec- 
tric and inside of it, respectively, in the case of the 
ideally plane dielectric surface with the equation of the 
profile -z = C (see eqs. (fT6|)- (|18|) \ The obtaining proce- 
dure for the potentials 1^1(2:), '/'2(z) and charge density 
n{z) is described in details in the previous section of 
the present paper, see eqs. ([26)) - ([59| . 

The next problem is concerned with the potentials 
Lpi[z,p) and ip2{z,p) obtaining. For these potentials 
one can use the eq. (|25p and the boundary condi- 
tions l|20p. (EH). In terms of the Fourier-transforms 
(fi{z,q} , (^2(2, q) over coordinate p of the potentials 
(piiz,p) and (p2{z,p) 



(,5i(z,p) = y d^gexp (iqp) (^i(z, q), 
(P2{z,p) = j (fqexp{i(ip) ^2(2, q) 



(60) 



the equations (f25|) have the following form: 

-'^i(^,q) 



d^ifiiz,q) 2 2dn{z) . 

- q ipi[z,q) = Aire - 



dz^ 



9^(^2(2, q) 
dz^ 



- g2(^2(z,q) = 0. 



(61) 



According to eqs. ([20]) . the boundary conditions 

concerned with these equations can be written as: 



{^i(z,q) - (^2(2, q)}^^^ = 

' d ip2{z) _ dipijz) 
dz dz 



e(q), 



-4n{eniz) + ^\ ^(q) - 4w(q; f) = ^^^^ 



2=« 



9(^i(z,q) d<f2{z,q) 
- e 



dz 



dz 



where C(q)j o'(q; ^^'^ the Fourier-transforms of the 
quantities ^(p) and (7(p;^), respectively (see eqs. (14), 
(18)): 

= / rf^9exp(iqp)|(q), 

(63) 

c^(p;0 = / d\exp(iqp)(5-(q;C). 

Let us consider the electric fields intensity perturba- 
tions caused by inhomogeneities of the dielectric sur- 
face rapidly decreasing at z ^ ±00. It is easy to see 
that the first equation in eq. (j6ip in general case can- 
not be solved analytically. But in two particular cases 
the analytical solution exists. In the first case, we solve 
eq. ([6T|) at z ~ ^ setting dn{z)/dfi equal to its value 
on the plane surface, z = ^: 



dn{z) dn{£_) 



(64) 



Such consideration is possible in the case when the typ- 
ical size of spatial inhomogeneities of the unperturbed 
charge density n{z) considerably larger than the typi- 
cal size of the spatial inhomogeneities of the potential 
1^1(2;, q) along z-axis: 



dn{z)\ d dn{z) 
dfi J dz dfi 



{^1(2, q)} 



-1 d(fiiz,(i) 



dz 



(65) 



Let us return to the discussion of the condition ([64)) 
below. 

Then, taking into account the assumption of rapidly 
fading field densities at z — )■ ±c», the solution of 
eq. ()60p can be given in the following form: 



ipi (z, q) = Ai (q) exp {~zb{q)) , 
(f2{z,q) = ^2(q)exp(zg) , 



where (see eq. ((64 




(66) 



(67) 



and Ai{q), A2{q) are obtained from the boundary con- 
ditions ([62)) . To this end, we put the expressions ([66) 



into the boundary conditions ([62]) and obtain the fol- 
lowing relations for the potentials (^i(z,q), (^2(2, q): 



(^i(z,q; 



exp{-{z-Ob{q)) 



sq + b{q) 
exp{(z~C)q) 



<^2(^,q) 
-An 



{[b{q) iE,{0-E2{0) 



eq + b{q) 

■ (eniO + {da{0/dO)] l(q; " 4w(q; } , 

(68) 

where (see eqs. ([22]), ((26)) ) 



( d(pi{2 



\ dz 



MO 



dz 



Let us emphasize that according to the boundary con- 
ditions (HOI), (HI]) (see also eqs. ([58]), ^) the values 
of the quantities -Bi(0; -^2(0 can be expressed as 



E^{0^Ea-2n\<j{0\>Q, 
E2{S.) = {E, + 2n\amie, 



(69) 



where the field intensity E^ is defined by the rela- 
tion ([M]) : 

Ea = E + iirens. 

Let us remind that the values of the potentials 
Lpi{z,q), (f2{z,q) at z — £^ do not coincide due to the 
fact that the potential continuity on the surface in the 
case of its inhomogeneous wavy structure is provided 
by the inequalities (see eqs. (fT6 |) -(|20 |) ): 



fiiz) I 



V2iz) I 



'^<y5i(^,q) = '^<^2(^,q), 



where (5.^1 (C, q) = ^{q) {d(piiz)/dz)^^^ + ipi{^,q). Ac- 
cording to eq. ((68l) . one can get: 



'^'^i(C,q) 



1 



[eqE2{0+b{q)Ei{0- 



eq + b{q) 

-in (en(0 + {d<7{0/dO)] l(q; ~ ^M^i; 

(70) 

where Ei{(^), i?2(C) ^-re still defined by the relations 
(69). It is easy to see from the obtained formulae ((68)) . 
([70]) that the gas of the volume charges can sufficiently 
affect on the potential of the electric field near the di- 
electric surface. 

Now let us show that the solution of eq. (pT|) in the 
forms (j66p . ((68)) is correct. As it is mentioned above, 
the condition of the existence of such solution is defined 
by the relation ((65|) . According to eq. ([67]) it can be 
expressed as follows: 



8(1 J dz 




The explicit expression for the derivative dn{^)/d^ 
can be obtained from eqs. ((35]), dM]), dU), (|46)) . 



(|52p , ([5^1 as in the case of degenerate charge gas above 
the dielectric surface, as in the case of nondegenerate 
one. In the second case the condition (|7T|) has a rather 
simple form: 

q^ > I3^e^ {{E + 2nensf + in^e^nl'j . (72) 

In the case, when the gas of charged Fermi-particles is 
degenerate at z <C zi (see eq. (ISS)) *) and low tempera- 
ture expansions (j52p take place, we can obtain the fol- 
lowing expressions for the volume charge density n{z) 
at ~ ^ and the electrochemical potential Tpo at z — 
(see eqs. dSH), dMD): 



.(Z) « 1^^3/2^ 



15 Ej - E^ 
32 nv 



2/5 



(73) 



where v and ip are still defined by the relations ([28|) 
with Q — —e, and Eq, E are expressed by (p6|) . By 
the use of the expression (|73p . one can write the con- 
dition ((7T|) : 

E^ - 52nensE - %n'^e^nl < 

, , I 15 (E -f 2nen,) \ (74) 
<C 4(7 e ' 

Accounting the solutions ((68)) are the Fourier- 
transforms of the potentials ipi{z,p), (p2{z,p) (see 
eq. ((SO)) ), the relations ((7T)) - ((7^ in general case are cor- 
rect for any value of g, including also the value q — Q. 
It is easy to see that the relation ((7^ does not satisfy 
such a requirement. The relation ((74p can take place 
at all values of q in the case of the external pressing 
field E that satisfies the following inequality: 



Q<E <E\ E' K 52nens 



(75) 



At E > E' the expressions (155)) do not take place. In 
this case, as in the case of the condition ((72)) realiza- 
tion, the equations ((6T)) must be solved by the use of 
the numerical methods. 

The case of the particular interest is the spatially pe- 
riodic inhomogeneities caused by the dielectric surface. 
As it is already mentioned in the present paper, such 
inhomogeneities are concerned with two-dimensional 
Wigner crystallization. In the most simple case of spa- 
tial periodic inhomogeneities the Fourier-transforms of 
the quantities £,{p), (y{p\0 (see eq. ((63))) can be ex- 
pressed in the form: 

^(q; = 2 51 ^(^"'^5 ^) + ^"'^) + (q - q^a)}, 

1 2 

C(q) = 2 51 ^^^<^^^ + +5{q- qa?)}, 

(76) 

where qaa {a — 1, 2) are the vectors of the recipro- 
cal two-dimensional lattice concerned with the spatial 



periodic charge distribution on the dielectric surface, 
Qaj (a = 1,2) are the vectors of the reciprocal two- 
dimensional lattice concerned with the spatial periodic 
wavy surface type, and d-{q^aa',£,) , CIQq?) ^'^e the am- 
plitudes of the corresponding surface heterogeneities. 
Of course, it is necessary to consider that the condi- 
tions (fT4|) take place, which in this case can be written 
as: 

qadiciai) < 1, gaal(qaa) < 1- (77) 



Then, putting the expressions ((76l) into eq. ((68|) and 
making inverse Fourier transformation according to 
eq. (|60|) . it is easy to obtain the following expressions 
for the potentials (^i(z,p), Lp2{z^p): 



-E^iO) + 4vr {eniO + da{0/dO k(qac) cos (q^^p) + 



exp(-(z-C)6((7Q^))^ 



^2 



exp((z-Ogag) 



b{qadiEliO- 



i (qajp) 



exp((z-^)q„^) ^ 
-47r > ^cr(qQ<T; cos (qa^p). 



eqaa + b{qaa) 



(78) 

The obtained expressions represent the solution of the 
potential distribution problem (so, the charge density 
distribution, too) in the area near the dielectric sur- 
face with the weak (see eq. (f77| ) spatially periodic in- 
homogeneities. Let us emphasize that for the expres- 
sions ([78)) validity it is no longer necessary to satisfy 
the conditions ([7T|) for all values of q. Its sufficiency is 
provided by the accomplishment of the conditions (|7T|) 
(or eqs. (|7^ . ([7i|) ) for two-dimensional reciprocal lat- 
tice distances qaa, qa^- 

Let us remind that we consider the simplest type of 
the spatial periodic inhomogeneities related to the di- 
electric surface. In the case of more complicated struc- 
ture of spatial-periodic homogeneities it is necessary to 
use the coefficients of two-dimensional Fourier series 
expansion for ^(p) and (t(p;^), which describe these 
inhomogeneities. 



5 Conclusion 



with the dielectric surface. The weak spatial inhomo- 
geneities caused as by the inhomogeneities of the sur- 
face itself, as by inhomogeneous charge distributions on 
it are taken into account. The case of "wavy" surface, 
in particular, the spatially periodic one, is concerned 
taking into account the possibility of the surface charge 
presence on it. The influence of the external pressing 
electric field acting on the system is also taken into 
consideration. It is shown that the presence of the gas 
of volume charges essentially influences on the value of 
the electric field potential in the area near the dielec- 
tric surface. Mostly, this fact plays an important role 
in the description of the deformation of the liquid di- 
electric surface caused by near-surface charges pressure 
on it. Authors of the present paper are working on this 
problem now. 

However, in our opinion, the solved problem is use- 
ful not only in two-dimensional Wigner crystallization 
aspect. The problem is worth concerning purely with 
the academic purposes as it can be related to the num- 
ber of classical problems of electrodynamics and sta- 
tistical physics. Due to this fact, in this paper we 
do not use the results of the real experiments on two- 
dimensional Wigner crystallization research. The for- 
mulations and obtained results in this paper can be 
used for the research of the influence of the volume 
charges near the liquid helium surface on the spatial 
inhomogeneous states of charges, which are adsorbed 
on the helium surface at the system parameters close 
to the experimental ones. 
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